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ABSTRACT 
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(N ' 

t— i \ Interferometers offer multiple methods for studying microlensing events and 

determining the properties of the lenses. We investigate the study of microlens- 
J> \ ing events with optical interferometers, focusing on narrow-angle astrometry, 

^ ; visibility, and closure phase. After introducing the basics of microlensing and 

interferometry, we derive expressions for the signals in each of these three chan- 

O ; nels. For various forecasts of the instrumental performance, we discuss which 

method provides the best means of measuring the lens angular Einstein radius 
#e, a prerequisite for determining the lens mass. If the upcoming generation 
Q-f of large-aperture, AO-corrected long baseline interferometers (e.g. VLTI, Keck, 

OHANA) perform as well as expected, 6e may be determined with signal-to-noise 
greater than 10 for all bright events. We estimate that roughly a dozen events 
per year will be sufficiciently bright and have long enough durations to allow the 
measurement of the lens mass and distance from the ground. We also consider 
the prospects for a VLTI survey of all bright lensing events using a Fisher matrix 
analysis, and find that even without individual masses, interesting constraints 
may be placed on the bulge mass function, although large numbers of events 
would be required. 



1. Introduction 

Gravitational lensing has been used for nearly a decade to study faint, compact masses 
in our galaxy. Although a large number of microlensing events have been detected, the lens 
masses and distances cannot (in most cases) be determined, meaning that only statistical 
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constraints on the lensing population may be derived from lensing surveys (e.g. Alcock et 
al. 1997, 2001). The determination of the lens mass and distance in individual events would 
be of great utility towards elucidating the nature of the microlenses. An example of this is 
the claim by Mao et al. (2002) and Bennett et al. (2002) that three long-duration events 
are likely massive black holes, with M ~ 10 — 3OM . Since individual masses could not be 
measured for these events, statistical arguments were employed to support the claim for large 
masses. If confirmed, Agol et al. (2002) have argued that these black holes would represent 
a new, significant population of black holes roaming the Galactic disk. 

The measurement of mass in microlensing events requires the determination of two 
quantities describing the event: (1) the lens-source relative parallax, 7Te, and (2) the angular 
Einstein radius, 6*e (Gould 2000). Measurement of the parallax, 7r E , requires the observation 
of the lensing event from viewpoints separated by > 1 AU. One way to do this is to observe 
the event simultaneously from the ground and from a satellite in solar orbit; the Space 
Interferometry Mission (SIM) is currently planned to do precisely this for a number of future 
microlensing events. It is also possible to measure 7Te for long-duration events, with event 
timescales £e ^few months, using the Earth's motion around the Sun to provide a distant 
vantage point. Since ~ 15 — 20% of lensing events towards the bulge have durations t = 2tg 
longer than 100 days, with high quality photometry the parallax may be measured for a 
significant fraction of events. 

The angular Einstein radius is extremely difficult to measure, since typical values are 
9e ~ O.5(M/O.3M ) 1,/2 milliarcsecond (hereafter mas), defying resolution by even the largest 
telescopes. In very rare cases (e.g. An et al. 2002), it is possible to measure #e during a 
caustic crossing, however generically 9e cannot be resolved by any single- aperture instrument. 
Resolution of sub-milliarcsecond scales typically requires long baseline interferometers, and 
future instruments like SIM should be able to determine 9 E for many events (Paczynski 1998; 
Boden et al. 1998a; Gould & Salim 1999) using narrow-angle astrometry. Unfortunately, 
SIM will not fly until later than 2009. Before SIM flies, a number of long baseline, highly 
sensitive interferometers will come online, for example the Keck Interferometer or the Very 
Large Telescope Interferometer (VLTI). 

These ground-based interferometers can study microlensing events in several different 
ways, which we will consider in this paper. Boden et al. (1998a) discussed how Keck and 
VLTI can measure 8e using narrow-angle astrometry. Delplancke et al. (2001) have pointed 
out that for massive microlenses, #e becomes comparable to the resolution \/B « 5 mas 
for B = 100m, A = 2.2/im. This allows the study of microlensing events via the partial 
resolution of the lensed images. One possible method, discussed by Delplancke et al. (2001), 
is the measurement of the decrement in fringe visibility as the microlensed images become 



- 3- 



resolved. In this paper, we also investigate the use of closure phase to determine the angular 
Einstein radius 9 E . Closure phase is free of many of the calibration issues afflicting visibility 
amplitude, however other concerns do arise. In the next section, we provide a review of the 
basics of microlensing. The following section gives an introduction to interferometry and 
closure phase. We then show how closure phase may be used to measure 6e, and compare 
our method to other proposed techniques. Since both Keck Interferometer and VLTI have 
already observed first fringes, this method promises to be an exciting technique for the 
determination of 9 E and thereby the lens mass, in the next few years. 



2. Review of microlensing 

In this section, we will focus on the case of a single lens and single (unresolved) source, 
since this describes the vast majority of lensing events. This situation was completely ana- 
lyzed by Paczyhski (1986), and we follow his notation here. Consider the geometry shown 
in Fig. 1. This shows a source at distance d s behind a lens at distance d t . The lens has an 
angular Einstein radius given by 



where M is the lens mass and di s = d s — di is the distance between lens and source. If the 
angular impact parameter between source and lens, in units of 9e, is u, then two images are 
produced on the sky along the lens-source axis at positions relative to the lens of 

x 1)2 = ^(«±v^+4). (2) 
These images are magnified relative to the intrinsic source brightness by factors 

A 12 = U ] + 2 ±1 (3) 

respectively. The total magnification is the sum of these two terms, 

A = A 1 + A 2 = =2A 1 -1 = 2A 2 + 1. (4) 

uyu 1 + 4 

In general, the lens and source exhibit finite relative proper motion, meaning that the 
impact parameter is a function of time. Usually, the relative motion is well approximated as 
linear in time, and can be written u(t) = w m in + A^Ae, measuring t with respect to the time 
of smallest impact parameter w m in- Here, £e is the event timescale, the time it takes the lens 



-4- 




Fig. 1. — (a) Illustration of the lensing geometry. Two images of the source are produced, 
at angular distances x\ and x<i from the observer- lens line of sight. The distances di, d s , and 
dig are between observer and lens, observer and source, and lens and source, respectively, 
(b) Illustration of the two images' motion vs. time. The dashed circle indicates the Einstein 
ring. The gray circles show the source at different times, while the light gray ellipses show 
the corresponding images. 



to move one angular Einstein radius relative to the observer-source line of sight, and is 
the unit vector on the sky along the relative velocity. Eqns. 2 and 3 then give £1,2 (£) and 
A li2 (t) as a function of time. An example is shown in Fig. 1. 

As noted by Boden et al. (1998a), the center-of- light of the two images shifts relative to 
the unlensed source position by 

A0coL = 8e^— (5) 
Also, note that at any given time, the separation between the two images is 

As = 6 E Vu 2 + 4, (6) 

and the ratio of magnifications is 

~A ± ~A + l~u 2 + 2 + uVu 2 + 4' 

These also become functions of time. 

It is apparent that the magnitude of the image motion is comparable to the angular 
Einstein radius. For a bulge event with M = O.3M , di = 4 kpc and d s = 8 kpc, this is 
roughly 9e = 0.55 mas, and scales like M 1 / 2 . In comparison, the resolution of the Hubble 
Space Telescope is ~ 100 mas, while the largest single aperture telescope (the 10m Keck) 
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has a resolution using adaptive optics of roughly ~ 50 mas. To have any hope of measuring 
#e, multi-element interfere metric arrays will be required. In the next sections we describe 
how interferometers may be applied towards microlensing. 

3. Review of Interferometry 

A stellar interferometer combines starlight from two or more apertures; the resulting 
interference fringes can be used to derive a great deal of information about the source being 
observed, including the source location (astrometry) and intensity distribution (imaging). 

3.1. Astrometry 

For observations with a finite bandwidth a fringe pattern will be formed when the 
optical pathlengths from the star, through the instrument, to the beam combination have 

— * 

been equalized. Given an aperture separation, called the baseline B, and a star in the 
direction s it follows from simple geometry that an additional delay d must be introduced 
into one of the arms of the interferometer, given by 

d = B-s + c (8) 

where c is any additional static delay internal to the interferometer. It follows that the 
differential delay Ad between two stars can be used to determine the angle As between 
them. Using laser metrology systems it is possible to measure internal delay differences 
at the nanometer level, which for a long baseline interferometer corresponds to astrometric 
precision at the micro- arcsecond level. This is the basis for interferometric astrometry, which 
has been demonstrated to achieve 100 micro-arcsecond precision (Shao et al. 2000), and is 
expected in the next several years to produce measurement precision approaching 10 micro- 
arcseconds (Boden et al. 1999; Delplancke et al. 2000). 

The practical aspects of very high precision interferometric astrometry are beyond the 
scope of this paper. However, the technique is quite challenging and requires a rather 
complex instrument design. The instruments currently being designed and built are not 
expected to begin operations until at least 2004-2005. In addition, due to the effects of 
atmospheric seeing, the very highest astrometric precision is likely only attainable by space- 
based interferometers. 
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3.2. Fringe Visibility 

The fringe pattern measured by an interferometer can be described as a complex quantity 
called the visibility (V\). It is related to the source brightness distribution on the sky (I\(s)) 
via 

fl(s)d(l [J> 

where A is the wavelength of observation and s points in the direction of the source. This 
relation is known as the van Cittert-Zernike theorem (Thompson et al. 1986; Lawson 2001). 
In effect, the source intensity distribution and the the fringe visibility measured by an inter- 
ferometer are related via a Fourier transform. 

In the case of a microlensing event, as long as the individual lensed images are unresolved 
by the interferometer, I can be modeled as two point sources located at s and si. 

I(s) = I S(s- s ) + hS(s- s\) (10) 

Defining the intensity ratio R — Ii/Iq and As = si — s*o we find 



V = TTH { 1 + H e"" A, " S ) (ID 



It is useful to write this visibility in a slightly different form 

V = \V\e i(t> (12) 

where 

\V\ 2 = (l + R 2 + 2Rcos(^As ■ B)^j (13) 

The quantity | | 2 (or simply V 2 ) is usually what is measured by an optical interferometer; 
it corresponds to the contrast of the observed fringes. The phase of the complex visibility is 
given by 

/ Rsm(^As-B) \ , 
(f) = arctan '— (14) 

VI + i?cos(f As-B)J 

However, the phase measured by a single pair of apertures is corrupted by the effects of 
atmospheric turbulence on very short timescales (in the optical and infrared ~ 1 radian in 
~ 10 milliseconds) and contains no useful information, unless it is measured simultaneously 
with a phase reference source - that is in fact identical to the astrometry approach discussed 
in the previous section. 

It is also important to note that V 2 measurements are subject to a wide variety of 
error sources, and that most of these errors introduce bias in the measurement (Tango & 
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Twiss 1980; Colavita 1999; Guyon 2002). That is to say that these error sources - which 
can include instrument vibrations, polarization mismatch, optical alignment errors, and the 
blurring effects of atmospheric turbulence - tend to reduce the measured fringe visibility 
amplitude. This makes calibration of measured visibilities quite challenging. Modern optical 
interferometers measure V 2 for one or a few baseline orientations; the precision attained is 
typically a few percent (Boden et al. 1998b) and in some cases as high as 0.3 % for bright 
stars under favorable conditions (Coude du Foresto et al. 2001). Unfortunately, it is not 
clear that the new generation of large-aperture interferometers equipped with AO systems 
and single-mode fiber spatial fibers will be able to attain such high visibility precision (Guyon 
2002). 



3.3. Closure Phase 

Despite the phase corruption introduced by the atmosphere, it is still possible to recover 
limited phase information without resorting to the technically complex phase referencing 
method, provided one interferometrically combines at least three apertures. In such a case 
one can form a quantity by multiplying the three complex visibilities formed over the three 
baselines. Since the visibility is the Fourier transform of the surface brightness distribution 
(c.f. Eqn 9), the product of three visibilities measured over a closed triangle is the surface 
brightness's bispectrum B(ki, fc 2 ), evaluated at &i = (ui, i>i) and k 2 = (u2, v 2 ) corresponding 
to the legs of the triangle. The phase of the measured bispectrum is known as the closure 
phase (Lawson 2001). The closure phase is immune to many forms of atmospheric corruption, 
which can be illustrated as follows: above each aperture there is a column of atmosphere 
with time-variable parcels of differing indices of refraction and hence optical pathlength. 
Thus the atmosphere above each aperture contributes a time-variable phase error, giving 

V = \v\ e i ^+^-<t>^ (15) 

where 4>i an d 02 are the phase errors associated with apertures 1 and 2 respectively, and 4>\i 
is an intrinsic phase associated with the source as measured by the 1-2 baseline. 

The bispectrum is thus 

t? 123 = I yd JV^I |^/ 3 | e i (^12+</ , 23+031+(</>l-</'2) + (02-</>3) + (03-</>l)) 

= |t4||^ 2 ||^3|e i(012+023+<?i3l) (16) 

We see that the atmospheric phase errors (as well as many other aperture-dependent phase 
errors) cancel. This is a well-known result, first applied in radio interferometry (Jennison 
1958). However, it is not immediately obvious what the closure phase represents. Below we 
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derive an expression relating the observed closure phase to the binary point source repre- 
senting a microlensing event. 

— * — * — * 

Assume 3 apertures, resulting in 3 baselines Bi,B 2 and B 3 . Note that 

B 1 + B 2 + B 3 = 0. (17) 

As before, we are looking at two point sources with intensity ratio R and separation As. On 
each baseline we measure a visibility V n given by Equation (9). 

V 123 = V X V 2 V 3 

= _L_ (i + ij e -¥^j (i + Re-^-B^ ^ + Re -^As-B^ 



(1 + i?) 3 

+jR 2 (e 2 F As.B 1 + e 2-A S -.B 2 + e ^As-B 3) + r3] (lg) 



the closure phase is thus 



(R 2 - R) Ei=i, 2 ,3 sin(^Ag ■ B L 
1 + R? + (R + R 2 ) E i=1)2 , 3 cos(f AS- B t 



123 = arctan [ - ^ - 7 _ ^ 7^ ^ I ( 19 ) 



There are a few things to note: the closure phase is always zero when R = 1, i.e. the 
source is symmetric. In addition, by Taylor expanding the sine terms and recalling Eqn. (17) 
it is easy to show that for the case when As 4 

As - 3 



^ X \JW\) (20) 

This implies that a source must be resolved by the interferometer in order to produce a non- 
zero closure phase; in the case of a partially resolved source the magnitude of the closure 
phase signal is very sensitive to the separation of the source components. 

Unlike V 2 , the phase measured by an optical interferometer is largely unbiased by mea- 
surement noise. In other words, the phase noise is zero-mean, and can be reduced by aver- 
aging over a sufficiently long time. The SNR for closure phase is given by (Tango & Twiss 
1980; Shao & Colavita 1992a) 



SNR, 



(iv^) +6 (^) Z + 4 (^) 



(21) 
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as compared to that for visibility and phase, 



NV 2 

sm ' = '' iTPV (22) 

In the photon-noise dominated regime (NV 2 3> 1), the SNR for closure phase and visibility 
scales as N 1 ^ 2 . However, for photon-starved sources, the SNR drops precipitously as N 3 ^ 2 , 
even worse than the oc iV scaling of the visibility SNR. Hence it is important to check whether 
microlensing sources will be photon-rich or photon-starved. 

In the if-band (2.2/im) a 15th magnitude source produces 0.16 photons cm -2 s -1 . Thus 
we would expect an 8-m class telescope to collect ~ 80,000 photons sec -1 . Assuming a 
coherence time r = 20msec, and a typical photon efficiency of 5%, we have N = 80 1 
per integration, safely in the photon rich regime when using modern low-noise detectors. In 
the absence of systematic errors a 3-element interferometer should achieve a closure phase 
precision of 10 -3 radians in 250 seconds. However, it remains to be seen if such levels can be 
achieved. To date, closure phase measurements have been made in the optical and near-IR 
by 2 groups. COAST achieves a closure phase precision of ~ 2 degrees (Baldwin et al. 1996), 
while NPOI achieves phase drifts of ~ 10 degrees hr -1 , which can be calibrated to the level 
of 1-4 degrees by looking at known single stars (Hummel et al. 1998). Recent progress in 
integrated optics - which due to their compact design are much less susceptible to systematic 
closure phase errors - may well yield a large improvement in the achievable accuracy (Berger 
et al. 2001). 

We suggest that a good method to correct for systematic closure phase errors would be 
to measure closure phase across a range of wavelengths. Given that many systematic error 
mechanisms give errors with a characteristic wavelength dependence that varies from that of 
the observed source, it should be possible to separate the signals. To illustrate, consider the 
case of a physical pathlength mismatch in the interferometric beam combiner, applied only 
to the light from one pair of apertures. Such a path mismatch will give a closure phase error 
that is inversely proportional to wavelength. Now consider the typical microlensing event; 
although the lens is achromatic, from equation (20) we see that the amplitude of the closure 
phase signal is inversely proportional to the cube of the wavelength. Hence, measuring the 
closure phase at several wavelengths should allow one to solve for both the source closure 
phase and the internal closure phase error. Clearly, this type of experiment would benefit 
from long baselines and wide bandwidth; the VLTI (200-m baseline) and AMBER (3-aperture 
combination from 1-2.4/xm) provide such a suitable combination (Glindemann et al. 2000). 
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4. Studying microlensing with interferometers 

Microlensing events are amenable to study by interferometers by all three methods de- 
scribed in the previous section. Microlensing causes an astrometric excursion of the center 
of light of the two images given by Eqn. (5), which may be measured using narrow angle 
astrometry (Paczyhski 1998; Boden et al. 1998a). Similarly, as the secondary image bright- 
ens, the two images become distinguishable from a single point source, causing a decrement 
in the visibility amplitude (Delplancke et al. 2001). Lastly, during the microlensing event 
the distribution of flux about the center of light becomes asymmetric, giving a signal in 
the closure phase. To measure 9e by astrometry, multiple measurements are required to 
map out the microlensing ellipse, while visibility and closure phase can give 9 E from a single 
measurement. 



4.1. Expected instrumental performance 

Now let us consider the measurement of these signals in practice. An important point 
to recall is that typical microlensed sources, by interferometrists' standards, are faint. This 
has several ramifications for the measurement of the signals in astrometry, visibility, and 
closure phase. As the application of these techniques to optical/IR wavelengths is still a 
relatively new field, there is considerable uncertainty regarding the measurement precision 
that will be achievable. Thus we will define a range of likely precisions (for the optimist and 
the pessimist) and evaluate the usefulness of the technique for each case. 

Astrometric measurements of faint sources require the presence of a bright (K < 13) 
phase reference within an isoplanatic patch, typically of order 20 arcseconds in the K band. 
The probability of a K < 13 star falling within 20" is roughly 17% (Shao & Colavita 1992b). 
Narrow-angle astrometry at the 100 micro-arcsecond level has already been demonstrated 
at PTI, and it is reasonable to expect the Keck and VLTI to produce the specified precision 
of 10-30 microarcseconds. However, narrow-angle interferometric astrometry is the least 
mature of the three techniques, and hence we will consider the utility of the technique for 
precisions of both 100 and 10 //as. 

Next consider measurement of the visibility amplitude for faint sources. The precision 
with which \V\ can be measured with small apertures (single- vq, where r$ is the Fried pa- 
rameter) can be as good as 0.3%. However, recall that in our estimate for the number of 
photons per coherence time, we assumed an aperture of 8m. Since a typical Fried parameter 
for K band is r ~ 60 cm, adaptive optics (AO) are required to correct distortions in the 
wavefront. Ideally, an AO system would completely flatten the wavefront, giving uniform 
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path delays across the aperture. This would restore the point spread function (PSF) to the 
familiar Airy pattern. Unfortunately, AO systems are not ideal, and in practice manage to 
get ~ 50% of the light into a tight core, and the remainder in a diffuse halo. This Strehl ratio 
sets the upper limit to the (squared) visibilities an interferometer can achieve, unless a spa- 
tial filter such as a single-mode optical fiber is used. However, such a spatial filter increases 
the maximum achievable visibility amplitude at the expense of losing light (max coupling 
efficiency ~ 78%, Shaklan & Roddier (1988)) and introducing visibility biases which depend 
on the Strehl ratio (Guyon 2002). In principle, these effects may be calibrated by observing 
point sources and monitoring the AO Strehl ratio in real time. In practice, such calibration 
is difficult and as yet unproven. At present, it is not known how well visibility amplitudes 
may be calibrated for large, AO-corrected apertures, but it may potentially be on the order 
of 5%. In summary, the expected precision in visibility amplitude (\V\) for the Keck and 
VLTI will likely be between 0.5 and 5%; we will consider the resulting measurement SNR 
for both cases. 

Lastly, we consider closure phase. As mentioned earlier, closure phase is free of many of 
the error sources afflicting visibility amplitude and phase. However, other systematics can 
arise which we have not anticipated. The current generation of interferometers have demon- 
strated closure phase precision at roughly the 1-degree level. However, there is no obvious 
reason why performance at the 0.001 radian (or 0.06 degree) could not be achieved, in partic- 
ular since the necessary high-precision phase- measurement equipment is being developed for 
differential-phase detections of hot exoplanets (Akeson, Swain & Colavita 2002). In addition, 
Segransan (2002) has discussed the possibility of detecting hot exoplanets using differential 
closure phase at the VLTI, an application which would require closure phase measurement 
precision on the order of 10~ 4 radians. Hence we will consider the closure phase SNR for 
both pessimistic (1 deg) and optimistic (0.001 rad) cases. 



4.2. Expected signal 

One of the differences between the study of microlensing via astrometry, as opposed to 
visibility or closure phase, is that the former does not rely upon resolving the two lensed 
images from each other. Because of this, the signal itself is independent of wavelength. 
However, a single astrometric measurement alone does not suffice to measure #e- in principle, 
at least 3 astrometric measurements are required to measure the microlensing signal and 
disentangle it from proper motion of the source. Since visibility and closure phase do involve 
resolving the two images, a single measurement does suffice to give #e, and the signal is a 
function of 9 E B/\. 
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Inserting Eqns. (6) and (7) into Eqn. (13) and expanding to lowest order in 9 E , we see 
that the visibility signal behaves as 



for #e <C X/2ixB. For a 3-element interferometer in an equilateral triangle, inserting Eqns. (6) 
and (7) into Eqn. (19) gives to lowest order 



again for 6>e <C X/ttB. Figures 2, 3 and 4 illustrate this. The steep dependence on the 
resolution means that there is a strong incentive to measure closure phase at the smallest 
possible wavelength. For example, the closure phase in J band (1/im) is more than 10 times 
larger than that in K band (2.2 /im). 

An actual microlensing event will produce time-varying signals in astrometry, visibility 
and closure phase, analogous to the familiar time-varying photometric lightcurves. Consider 
a typical microlensing event with a O.5M lens at distance d\ — 6 kpc, giving 9e = 0.43 
mas. Assume u min = 0.35, so that the maximum magnification is A = 3, and the maximal 
brightness ratio of the images is R = 1/2. An interferometer with a 100m baseline operating 
at 2.2 fim (fT-band) has a resolution of 4.5 mas. If we have three such baselines arranged in 
an equilateral triangle, they would observe the signals in astrometry, visibility, and closure 
phase shown in Fig. 3. The astrometric signal has a maximal amplitude of ~ 0.35#E) which 
in this example is about 0.15 mas. The visibility, which is usually minimized at closest 
approach (u = Mmm), is here |V| = 0.85. The closure phase has a maximum value of 
0123 = 1-5°. With present technology, precisions of ~ 5% in |V| and ~ 1° in closure phase 
may be achieved, so at present such events should be studied by visibility. However, advances 
in measurement precision and observation at shorter wavelengths can significantly enhance 
the signal-to- noise. In the next subsection, we consider a reasonable range for the expected 
instrumental performance and show how this affects the prospects for determination of #e 
using interferometers. 




(23) 




(24) 



4.3. Prospects for mass measurement 



Except for highly resolved events, the astrometric signal scales as ~ E , while the 
visibility signal scales as ~ Q\ and the closure phase as ~ Q\. These steep scalings improve 
the SNR in the derived Einstein radius; for example a SNR of 5 in visibility translates into a 
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Fig. 2. — The first panel shows the visibility amplitude as a function of u for varying degrees 
of resolution: 9eB/\ = 0.05, 0.07, 0.1, 0.14, 0.2 from top to bottom. The second panel shows 
the corresponding closure phase, in degrees. Note that for highly resolved events, the signals 
do not follow the approximate Eqns. (23) and (24). 




Fig. 3. — Interferometric signal for a microlensing event with 6*e = 0.43 mas and u m i n = 0.35. 
The first panel shows the magnitude of the astrometric motion of the center of light. The 
second panel shows the visibility amplitude, and the third panel plots the closure phase. We 
have assumed a three element interferometer in an equilateral triangle, with arms of 100 m 
each, observing at 2.2 /im. We also assume that Earth rotation allows a maximal projection 
of the image separation along at least one of the baselines. 
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SNR of 10 in E , while a closure phase SNR of 5 becomes a #e SNR of 15. Figure 4 plots the 
SNR in #e from astrometry, visibility and closure phase, assuming respectively optimistic 
errors (left panel) of 10pm in astrometry, 0.5% in \V\, and 0.001 rad in closure phase, and 
(right panel) pessimistic errors of 100 pm astrometry, 5% visibility amplitude, and 1° closure 
phase. If errors of 0.001 can be achieved in J band, then events with #e > 0.1 mas can be 
measured with SNR> 10. 

The next step is to determine the brightnesses and Einstein radii for typical events. 
Figure 5 illustrates the properties of bulge events. The first panel shows the distribution 
of peak / magnitudes (or in cases where the peak was not observed, the brightest observed 
magnitude) for events observed by the OGLE collaboration during the 1997-1999 seasons. 
Assuming that sources with K < 14 can be observed interferometrically, and assuming that 
I — K — 1.5, we estimate that roughly 20% of events are bright enough to be followed up 
with interferometers. However, future surveys may go deeper than OGLE-II, returning a 
smaller fraction of bright events. Next, we turn to the distribution of Einstein radii. Since 
this has not been measured, we must estimate theoretically what the distribution will be. 

The rate of microlensing events depends upon the properties of the lenses and source 
stars. We assume two types of stars, bulge stars and disk stars. For the disk stars, we use 
the mass function of Gould et al. (1997), namely dN/dM oc (M/M b ) a , with M b = O.59M , 
and a = -0.56 for M min < M < M b , and a = -2.21 for M b < M < M max . We take 
M m i n = O.1M and M max = 100M Q . The disk density profile is taken to be p(r{) = poe ri ^ rd , 
where po = 8 lO 7 M /kpc 3 is the local density, = 3.5 kpc is the disk scale length, and 
r\ again is the distance of the lens from us towards the source, taken to be at the Galactic 
center. The velocity distribution of disk stars is taken to be a flat rotation curve v c = 220 
km/s, along with velocity dispersion in each transverse direction of a d = 30 km/s. For 
the bulge stars, we use the mass function of Zoccali et al. (2000), which has the same 
form as the disk MF but with M b = 1M , a = -1.33 for O.15M < M < M b , and 
a = —2 for M b < M < 100M Q . The bulge density distribution we use is based on the 
barred model of Han & Gould (1995), with p(ri) = poexp[— ^ (r s — r/) 2 ], with central density 
Po = 2.07 10 9 M Q /kpc 3 , w = (cosi/xo) 2 + (sini/i/o) 2 with scale lengths xq = 1.58 kpc, 
yo = 0.62 kpc and inclination % = 20°. Bulge stars are taken to have no net rotation and a 
velocity dispersion in each transverse direction of <r b = 110 km/s. 

Lenses are drawn from both the disk and bulge populations. For simplicity, we assume 
that all sources lie at the galactic center, with the bulge velocity distribution. We neglect 
the motion of the local standard of rest relative to the flat rotation at v c . The distribution 
of Einstein radii is given in the Appendix by Eqn. (A8) and Eqn. (A7). Using the above 
parameter values, the resulting rate distribution is plotted in the second panel of Fig. 5. 
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Fig. 4. — SNR in E . The solid line corresponds to the signal-to-noise ratio in 6 E derived 
from closure phase, the dashed line corresponds to visibility, and the dotted line corresponds 
to astrometry. The left panel shows 'optimistic' signal-to-noise ratios, assuming 10 //as errors 
in the astrometry, 0.5% errors in |V| and 0.001 rad errors in the closure phase. The right 
panel shows 'pessimistic' signal to noise using 100 /xas astrometric errors, 5% errors in |V| 
and 1° errors in closure phase. We assume astrometric measurements are taken at u — 1.4, 
closure phase is measured at u = 0.86 in J band (1 jum) and visibility is measured at u — 0.35 
in K band (2.2 /jm) where higher Strehl ratios may be achieved. 




100 



Fig. 5. — The first panel shows a histogram of / magnitudes for 213 microlensing events 
observed by OGLE during 1997-1999. Of these, 26 have I < 15, and 48 have I < 15.5. 
The second panel shows the distribution of lensing rate with respect to angular Einstein 
radius, dT / d\ogO E . The solid line shows the total rate distribution, the dashed line shows 
the distribution of events with t E > 50 days, and the dotted line shows that for events with 
t E > 100 days. 
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Clearly, a large fraction of events should have 8e accessible to instruments like the Keck 
Interferometer or VLTI. 

As mentioned earlier, for very long timescale events (t E > 50 days) the parallax may be 
measured with sufficiently high quality photometry, meaning that measurement of 9e gives 
the mass and distance. Estimating that ~ 15% of events are of sufficiently long duration, 
that ~ 10% are bright enough to be observed with VLTI or Keck, and that 1000 events 
are detected every year we expect to measure masses for about 15 microlensing events from 
the ground every year. This should be interesting for constraining the high end of the mass 
function. 

However, VLTI and Keck can observe many more events than just those showing par- 
allax, if there is justification to do so - for example, if such measurements could better 
constrain the lens properties. Let us wildly speculate about a massive campaign from VLTI 
to observe 100 events per year. The basic observable from such a campaign is the joint distri- 
bution of event timescales £e and angular Einstein radii 9e- An expression for dT /g^eg^e is 
derived in the appendix in Eqn. (A10). The resulting rate distribution is plotted in Fig. (6). 

A massive survey can measure a noisy sample of this distribution. To estimate how 
well these measurements may be used to constrain the underlying lens properties, we use 
the Fisher matrix (Gould 1995; Tegmark et al. 1997). For a likelihood distribution L (here 
dT /g^eg^e) with parameters Aj, the covariance matrix of the parameters estimated by sam- 
pling L is given by 

Cij = (h'^j (25) 

where the matrix b has elements 2 

/ dlogL<91ogL \ 

h > = \ dx t dx, ) (26) 

and the angular brackets denote integration over possible observables (here 6e and tE) 
weighted by L. We vary eight parameters: the bulge velocity dispersion c% mass func- 
tion parameters a±, a<i and m p , and the same parameters for the disk lenses. Results are 
listed in Table 1. We find that the corresponding disk and bulge parameters are somewhat 
degenerate, as may be expected, but that large numbers of events can allow constraints to be 
placed on the bulge parameters. For example, with N = 400 lenses the error on the slope of 
the high end of the bulge mass function becomes Aa 2 = 0.2, and the error on m p b becomes 
O.18M , sufficient to detect a break in the mass function and to localize the break position. 



2 We follow Gould (1995)'s notation for the Fisher matrix, which differs from Tegmark ct al. (1997)'s 
notation by a partial integration. 
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The caveat to this statement is that we have made several simplifying approximations, for 
example not exploring degeneracies between the mass function parameters and the density 
profile parameters. 

5. Discussion 

Interferometry is poised to revolutionize the study of microlensing events. Until now, 
microlensing has suffered the difficulty that masses of individual lenses cannot be measured, 
severely limiting the information able to be extracted from lensing surveys. The problem 
has been that the two quantities needed to measure mass and distance, the relative parallax 
7Te and the angular Einstein radius 9 E , are not regularly measured. The parallax may be 
measured for long duration events with high quality photometry, however measurement of 
#e requires resolution on the order of a milliarcsecond, necessitating interferometers. The 
upcoming Space Interferometry Mission (SIM) should measure masses and distances for a 
large sample of lenses, answering the question of the microlenses' nature. Well in advance of 
SIM, however, ground-based interferometers can also provide useful measurements of lensing 
events. 

As mentioned earlier, microlensed sources are generally much fainter than the typical 
sources studied by optical interferometers, meaning that large apertures (~ 8m) are required. 
Both the Keck Interferometer and the VLTI can measure visibility amplitude for microlensing 
events using their largest apertures, but only VLTI can measure closure phase using 3 large 
apertures; the Keck would be required to employ one of the 1.8 m outrigger telescopes which 
collect considerably fewer photons. The signal measured by interferometers is a function of 
the Einstein radius in units of the resolution, 9-&B/X. Because of this, there is great advantage 
to go to shorter wavelengths. However, shorter wavelengths require the use of adaptive optics 
(AO) systems. Since AO makes the accurate calibration of visibility very difficult, but has 
a smaller impact on the calibration of closure phase, there are obvious advantages to using 
closure phase. For a numerical example, a microlensing event with 9 E = 0.5 mas could be 
observed at 10 /im without AO, but 9 E B/\ ~ 0.02 giving a visibility signal of V xs 0.99, 
which would be extremely challenging to distinguish from a point source. The same event 
observed at 2.2 /im using AO has 9 E B/\ ~ 0.11 giving a closure phase signal 0123 ~ 2.5° 
which can already be measured. 

Only a small fraction of events are expected to be bright enough (K < 14) to be observed 
interferometrically. However, certain fainter events may also be accessible to interferometers. 
If a bright (K < 13) star falls within the isoplanatic angle, then phase referencing may be 
employed to extend the coherence time significantly. As noted earlier, for sites with small 
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Fig. 6. — Contours of cir/iilog^E^logtE! spaced by decades. 
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value 


error (N = 




110 km/s 


108 km/s 


m Pjb 


1 M 


3.6 M 




-1.33 


7.9 




-2 


3.8 




30 km/s 


466 km/s 


m p4 


0.59 M 


4.0 M 


aid 


-0.56 


43.3 


a 2 d 


-2.21 


9.3 



1) 



Table 1: Errors on parameters from measuring dT / dOvdtv- The second column lists the 
adopted values for each parameter, and the third column lists the measurement errors esti- 
mated using the Fisher matrix. The errors values are to be multiplied by N^ 1 ^ 2 , where N is 
the number of bright lens events observed by the interferometer. 
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isoplanatic patches the probability of finding a suitable bright star is poor. Additionally, 
this technique is quite complex and as yet unproven, however in principle this could allow 
the study of microlensing events as faint as 20th magnitude, reaching the bright end of LMC 
events. Phase referencing must be employed to perform narrow angle astrometry; our results 
indicate that events for which phase referencing is possible may be more profitably studied 
with visibility or closure phase. 

We expect that ~ 15 events every year will be bright enough (K < 14) and have suffi- 
ciently long duration (tg > 50 days) to permit the measurement of mass and distance. We 
have shown that a fairly large fraction of events accessible to ground-based interferometers 
should allow measurement of #e with high signal to noise. We also investigated the prospects 
for a massive follow-up campaign by VLTI, and found that statistical information on the 9e 
distribution, even without individual mass measurements, can allow constraints to be placed 
on lens properties like their mass function. 

Even if our estimates turn out to be overly optimistic, interferometers will still be able 
to elucidate the nature of claimed black hole candidates (e.g. Mao et al. 2002; Bennett et 
al. 2002). Agol et al. (2002) have suggested that current microlensing data indicate the 
presence of a significant population of intermediate mass black holes roaming the Galactic 
disk; interferometers will be able to confirm or reject this possibility. 

In this paper, we have focused on ground-based interferometers, however our results 
apply also to space-based interferometers like the Space Interferometry Mission (SIM). SIM 
is primarily an astrometric instrument, however it can also measure fringe visibilities. Nomi- 
nally, the target precision expected for SIM is 1% in V 2 (M. Shao 2002, priv. comm.). SIM's 
baseline is 10m, and typical wavelengths are A ~ 0.6/im, giving resolution of about 12 mas. 
Hence, SIM can determine #e with SNR> 10 from visibility alone, entirely independently of 
the astrometric determination, for events with #e > 0.44 mas. From Figure 5 we see that 
this comprises a large fraction of the events. The visibility measurements come for free with 
the astrometric measurements, and should significantly increase the precision of SIM mass 
measurements, as long as effects such as crowding do not pose too great an obstacle. In 
addition to measuring # E , SIM also determines 7Te, the lens parallax, by measuring the time 
of the peak of the photometric lightcurve. Since the peak of the visibility signal coincides 
with the peak of the photometric signal, SIM visibility measurements could also be used to 
determine 7Te. However, since the variation in 1 — V 2 is so shallow near the peak this method 
may not prove to be as precise as ordinary photometric parallax measurement. 

One of the most (potentially) exciting prospects is a topic we have not discussed in this 
paper, binary microlensing. For a binary lens system, complicated caustic structures can 
arise, leading in favorable cases to the production of 5 images of the source. For a spectacular 
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example of this, see An et al. (2002). During caustic crossings, the magnification can get 
exceptionally large, e.g. factors of 30, making these events bright enough to observe with 
interferometers. The five images are currently unresolvable from each other, however VLTI 
and Keck offer the prospect of imaging the multi-image pattern. With a multi-aperture 
system (required for closure phase), one obtains several visibility measurements and one 
closure phase at the same time, possibly allowing the reconstruction of complex events such 
as caustic crossings. 
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A. Appendix - Microlensing rate 

The microlensing optical depth r = f nadl is given by 
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where dN/dM is the mass function of lenses, normalized by J(dN/dM)dM = 1, n(ri) = 
p(ri)/M is the number density of lenses, M = J (dN/dM)MdM is the average mass and p 
is the mass density of lenses. From Eqn. 1, we then have 
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The optical depth distribution describes the instantaneous properties of lensing events 
at any given time, but we are interested in the distribution of all events. This is described 
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by the lensing rate T = ncrv. If the transverse velocity distribution of source stars is f s (v s ), 

and that of the lenses is fi(vi), we can write the lensing rate as (Griest 1991) 

r (IN r r ° 
V ----- I dM— j <lrf2n() v nin)rr{n) 



dM 



(A3) 



where 



v T = Vl - (1 - x)V - xv s . 



(A4) 
(A5) 

Here, v Q is the observer's velocity transverse to the line of sight, and x = r t /r s . We assume 
all sources are bulge stars, with no net rotation and a velocity dispersion in each transverse 
direction of cr b . Lenses are assumed to reside either in the bulge or in the disk, with the 
latter population possessing a flat rotation curve v c , along with velocity dispersion in each 
transverse direction of a d . For bulge lenses, this gives 
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As usual, the angular integral gives a Bessel function, but note that it is a modified Bessel 
function due to the lack of an i in the exponent. Thus we have 
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For disk lenses, the expression is similar, with (1 — x)v c —r juu c &i±u. v ± -r ^ ju b 

From Eqn. (A3), the distribution of rate with respect to angular Einstein radius is 
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Lastly, we also want the joint distribution of rate with respect to 9 E and t E . This is 



dT 



d9 E dt E 



pr 

Jo 



dN 

dririn(ri)m(9 E ,ri)—[m(9 E ,ri)] \ d 2 Vid 2 v s fi(vi)f s (v s )v T 6[t E - 



t E 



(A9) 



-22 - 



For bulge lenses, this is 
dT 



d9 E dt E 
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For disk lenses, the expression is similar, with (1 — x)v c — > xv c and (1 + x 2 )a 2 



a 2 + x 2 a 2 . 
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